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We compute low-lying eigenmodes of the gauge covariant Laplace operator on the lattice at finite 
temperature. For classical configurations we show how the lowest mode localizes the monopole 
constituents inside calorons and that it hops upon changing the boundary conditions. The latter 
effect we observe for thermalized backgrounds, too, analogously to what is known for fermion 
zero modes. 

We propose a new filter for equilibrium configurations which provides link variables as a truncated 
sum involving the Laplacian modes. This method not only reproduces classical structures, but also 
preserves the confining potential, even when only a few modes are used. 
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1. Introduction 

We study low-lying eigenmodes of the gauge covariant Laplace operator 

4 

-A(j) n (x) = KU*) , A# = £ <(x)5,+a,v + Ul a \y)8 x ^ y - 28 ab 8 x 



(1.1) 



with Ufi(x) a given lattice configuration in the fundamental (or adjoint) representation of SU (2). We 
will use these modes as an analyzing tool, like it has been done with fermionic (near) zero modes. 
The latter owe their existence to the index theorem, and in smooth cases they are localised to 
the nontrivial background. Laplacian modes are computationally cheaper because they are space- 
time scalars and do not suffer from chirality problems nor doubler problems. They seem not to 
be directly related to topology, but it has been observed that they are sensitive to the location of 
instantons [|l|, ||]. The scaling properties of their localization have been investigated recently |gp. 

Here we will present two ideas concerning Laplacian modes [Q]. The first one is to study what 
can be learnt from the profile of the modulus of the lowest mode (the 'ground state probability 
density') \(f)o(x) j 2 , in particular with changing boundary conditions. The second idea is to introduce 
a novel low-pass filter. It concerns the reconstruction of the gauge background from a few Laplacian 
modes. The purpose of both methods is to identify the underlying infrared degrees of freedom in 
the gauge field that are responsible for features like confinement. 

2. Profile of the ground state density 

2.1 Caloron backgrounds 

As a testing ground for the Laplacian modes we first investigate a caloron of maximally non- 
trivial holonomy with its monopole constituents [Q, ^] put on a 16 3 -4 lattice 1 . Fig. [I] (left) shows 
the action density along the line connecting these monopoles, which are clearly visible as two self- 
dual (and almost perfectly static) lumps. As another signal the Polyakov loop goes through 1 2 and 
— 1 2 - which amounts to a local symmetry restoration - at the monopole cores. 





Figure 1: A caloron with its constituent monopoles. Left: the action and topological density (not dis- 
tinguishable, both multiplied by a factor 400) and the Polyakov loop. Right: the modulus of the lowest 
Laplacian mode with periodic (boldface line) and antiperiodic (thin line) boundary conditions. 



This is done by calculating links from the continuum gauge field, followed by a few steps of cooling. 
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The ground state of the Laplacian in this background is shown in Fig. [I] (right). One can see 
that the presence of monopoles is reflected by a maximum resp. a minimum in the profile of the 
lowest mode (although shifted by up to 2 lattice spacings). In addition to the mode periodic in 
Euclidean time we have depicted the antiperiodic one. In this mode the role of the monopoles are 
interchanged, which can be understood from a symmetry of the caloron. Hence the lowest-lying 
Laplacian mode 'hops' as the result of changing the boundary conditions. It behaves similar to 
the caloron fermion zero mode [{7]] in the context of which complex boundary conditions were 
introduced first. 

2.2 Thermalized backgrounds 

In this subsection we explore the Laplacian modes as an analyzing tool on thermalized gauge 
field configurations. Fig. ^| shows the lowest-lying mode in a generic background obtained on a 
16 3 -4 lattice at j8 = 2.2 (which amounts to T = 0.75 T c ) confirming the expectation that these 
modes are free of UV fluctuations. We allow for a complex phase in the boundary condition, 
0(^4+^) = exp(27n'£)0(jt4) parametrised by an angle £ G [0,1/2] (the remaining £'s can be 
mapped to this interval by charge conjugation). The 'evolution' of the Laplacian mode with this 
angle is shown vertically in Fig. ||, while the 3 rows represent different but fixed lattice planes, 
which contain the global maximum of the mode. In the example we find 3 lattice locations which 
become the global maximum within some £ -interval. Therefore, the lowest Laplacian mode in a 
thermalized background is hopping, too. 

At the hopping points in £ the value of the maximum 2 as well as the inverse participation ratio 




Figure 2: The lowest Laplacian mode in a thermalized background with £ = 0.02,0.355,0.48 from left to 
right and in lattice planes (x2,X4) = (12, 1), (4,4), (6,3) from top to bottom. The vertical scale is 0.05. 



The global minimum is not stable enough to employ it for analyzing purposes. 
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(taking values between IPR = 4 and IPR =11) decrease and the mode rearranges itself. For some 
boundary conditions the lowest Laplacian mode can be characterized as a global structure. 

Analyzing 50 independent configurations the lowest mode hops up to four times as a function 
of £. Apart from the Laplacian mode being wave-like, the hopping phenomenon resembles the 
behaviour of the fermion zero mode However, in both cases it is not obvious which gluonic 
feature, say in the topological charge, discriminates the preferred locations. Interestingly, we have 
found a correlation of the maximum of the periodic and antiperiodic Laplacian mode to positive 
and negative Polyakov loops, respectively. This is the same tendency as for calorons (see Fig. [[J) 
which could be understood as an argument for calorons 'underlying' the quantum configuration. 
Alternatively, the Polyakov loop could play a role defining pinning centers for the Laplacian mode 
in the spirit of Anderson localisation in a random potential 

The lowest Laplacian mode in the adjoint representation we have found to have minima at 
the caloron monopoles (see also [|l|, [10|]) that extend to two-dimensional sheets for antiperiodic 
boundary conditions. In thermalized backgrounds, the maxima of the adjoint modes are correlated 
to the fundamental ones with the same boundary condition, but are stronger localized than the latter. 



3. A Fourier-like filter 



Laplacian ('harmonic') eigenmodes can be used to define a new Fourier-like filter which we 
present now [fflj. We were inspired by the representation of the field strength in terms of fermionic 



modes in [11], but here we aim to reconstruct directly the link variables on which any observable 



can be measured. To this end we combine the definition of the lattice Laplacian, Eq. (|L1|), with a 
spectral decomposition and at y = x + jx immediately obtain 

U?(x) = - £ ^<(jc)0„* fc (jc + A) , /=««'2 (3.1) 

n=\ 

The idea is to truncate the sum on the right hand side of this equation at some N <C JV. The question 
arises, how to relate such an expression to a unitary link variable. Here the charge conjugation 
helps, since it guarantees that every eigenvalue is two-fold degenerate with eigenfunction related 
as <t>' a (x) = e ab ty* b {x). It follows that the corresponding bilinear expressions in Eq. (3.1) add up to 
an element of SU (2) up to a factor. The same applies to the staple average in cooling or smearing. 
We divide by the square root of this factor, which is the determinant and positive for all practical 
purposes, and obtain the final filter formula 

U« b (x) N = (- £ K [^{x)^\x+^)-e^ n {x)^{x + ^)e db ^ /Vd^TO (3.2) 

The quality of the filter is controlled by N, where N = JV reproduces the original configuration 
exactly. In the other extreme case, iV= 1, the filtered links U^(x) can be shown to be pure gauge. 
From the behaviour of <p a (x) under gauge transformations it follows immediately that U^x) trans- 
forms as a link (and no gauge fixing is involved in the filter). 

The filtering can be performed with Laplacian modes of any boundary condition. Then Eq. 



(3.2) is just the average over opposite boundary conditions £ and — £. This additional parameter 
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Figure 3: Action and topological density and Polyakov loop of the caloron of Fig. [T] (left) filtered with N = 4 
(left, densities multiplied by 100) and N = 150 (right, densities multiplied by 400) modes (£ = 1 /4). 

completely fixes the Polyakov loop at N = 1 to tr &(x)/2 = cos (2nQ, while for nontrivial cases 
N >2 the Polyakov loop is observed to fluctuate around this value. It follows that in order to 
optimize the filter for our circumstances (confined phase, nontrivial holonomy calorons) the best 
choice is £ = 1/4, i.e. halfway between periodic and antiperiodic boundary conditions. 

In order to test the filter we again start with the caloron background. Fig. ||] shows the action 
and topological density as well as the Polyakov loop measured on the filtered links U^(x), to be 
compared with Fig. |l| (left). The filter with the number of modes as low asN = 4 starts to reproduce 
the classical structures qualitatively, while for N = 150 the agreement is almost perfect. 

What is even more remarkable is that the filter method preserves the string tension. We plot 
the logarithm of the Polyakov loop correlator measured on 50 configurations at T = 0.75 T c in Fig. 
^ It reveals a clearly linear behaviour with a slope reproducing the original one within 15%. The 
minimal number of modes, N = 2, suffices for this behaviour. We conclude that the confining 
properties of lattice gauge theory are captured by the lowest Laplacian modes. The Polyakov loop 
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Figure 4: Polyakov loop correlator plotted on an inverse logarithmic scale over the distance for different N 
(at £ = 1/4) compared to the original one. 
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correlator after filtering has no sign of a Coulomb regime, since the filter has washed out short 
range fluctuations. At zero distance the values ((tr &(x)) 2 ) are smaller than the original one which 
indicates that the distribution of Polyakov loops over the lattice sites has become narrower. 

Since the filter is not biased to classical solutions nor to particular degrees of freedom in 
the gauge field (like monopoles or vortices), it is interesting to have a closer look at the vacuum 
structures that emerge when the filter is applied to generic configurations. In both the Polyakov loop 
and the action density, more and more fluctuations appear with an increasing number of modes kept 
in the filter. It might well be that the fluctuation at N = 2 are the minimally required ones for the 
long-range physics. 

Concerning the action density, these objects are isolated peaks (which are non-static and not 
necessarily (anti)selfdual). This phenomenon is counterintuitive as the filter uses the lowest-lying 
Laplacian modes, however, this is also the case for the classical calorons, see Fig. [| (left). We have 
found a correlation of the peak structure of the filtered action density to that emerging after cooling 
or smearing in an early stage. More work has to be done to better understand the (A^-dependent) 
spiky structures induced by the filter. 
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